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Abstract 

The energy spectrum of cosmic-ray antiprotons (p's) has been measured by BESS 
successively in 1993, 1995, 1997 and 1998. In total, 848 p's were clearly identified 
in energy range 0.18 to 4.20 GeV. From these successive measurements of the p 
spectrum at various solar activity, we discuss about the effect of the solar modulation 
and the origin of cosmic-ray p's. Measured p/p ratios were nearly identical during 
this perfiod, and were consistent with a prediction taking the charge dependent 
solar modulation into account. 
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1 Introduction 



The origin of cosmic-ray antiprotons (p's) is still holding a fundamental ques- 
tion even after their spectrum has been measured precisely in our previous 
work [1]. In recent solar-minimum period, we have collected 458 p's in the en- 
ergy range between 0.18 and 3.56 GeV and observed a clear peak in the spec- 
trum around 2 GeV, a generic feature of "secondary" p's which are produced 
by the interaction of Galactic high-energy cosmic-rays with the interstellar 
medium. It has become evident that the dominant component of cosmic-ray 
p's is the secondary p's. While the energy spectrum and absolute flux are well 
reproduced by several theoretical calculations [2,3], there remains some diver- 
sity among those calculations in the low energy region. They seem to arise 
from the different treatment of tertiary interactions of p's with the interstellar 
medium, p production mechanism such as nuclear sub-threshold effects, galac- 
tic propagation and solar modulations. Moreover, although the data have still 
large statistical errors at low energies, the observed flux appears to be higher 
than those calculations in the low energy region below 1 GeV where the sec- 
ondary flux decreases sharply. Therefore, we can not rule out the primary 
origins, such as annihilation of neutralinos [2,4,5] or evaporation of primordial 
black holes [6,7]. 

The cosmic-ray particles are decelerated in entering the heliosphere by the 
expanding solar magnetic disturbances, and their energy spectra observed at 
1 A.U. should be deformed. The flux variation is expected to be small for 
secondary p's at low energies because of their relatively hard spectrum, while 
any softer component would exhibits a large variation. It has been pointed 
out that additional soft component from a novel source can be singled out 
by measuring the p spectra at different solar periods [5]. It has been recently 
suggested that the solar magnetic polarity and the charge dependent effect 
cannot be neglected and the modulation follows a 22-year cycle rather than 
an 11-year cycle [3]. Actually, measurements of e~/He 2+ ratio [8] and e~ /p + 
ratio [9] variations over a long period appear to support this charge depen- 
dence. Since protons and antiprotons are different only in their charge sign, the 
simultaneous measurements of both particles are most suitable for studying 
the charge dependence of solar modulation. 

We report here a new measurement of cosmic-ray p spectrum in the energy 
range between 0.18 and 4.20 GeV, based on 384 events detected by the BESS 
spectrometer in 1998 flight. In total, 848 p events has been accumulated from 
the measurement in successful four flights in 1993, 1995, 1997 and 1998. The 
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spectrum measured in 1998 is compared with the previous results [1,10-14], 
in order to observe changes with solar activity. 



2 BESS Detector 



The BESS detector, shown in Fig.l, was proposed and developed as a high- 
resolution spectrometer to perform searches for rare cosmic-rays as well as 
precision measurements on absolute fluxes of various cosmic-ray particles [15- 
17]. Various new detector technologies developed for collider experiments were 
incorporated to the spectrometer. A uniform field of 1 Tesla was produced by 
a thin (4 g/cm 2 ) superconducting coil [18], through which particles could pass 
with small interaction probability. The magnetic-field region was filled with 
central tracking devices. This configuration achieved a geometrical acceptance 
of 0.3 m 2 Sr, and it was an order of magnitude larger than those of previous 
cosmic-ray spectrometers for balloon experiments. The rip tracking in the cen- 
tral region was performed by fitting up to 28 hit-points, each with 200 jum 
resolution, in a jet-type drift (JET) chamber and two inner drift chambers 
(IDC's) resulting in a magnetic-rigidity (R = Pc/Ze) resolution of 0.5 % at 
1 GV/c. Tracking in the z coordinate was done by fitting points in IDC's 
measured by vernier pads with an accuracy of 470 /xm and points in the JET 
chamber measured by charge division with an accuracy of 2.5 cm. The con- 
tinuous and redundant 3-dimensional tracking enabled us to recognize multi- 
track events and tracks having interactions or scatterings, thus minimizing 
such backgrounds. The upper and lower scintillator hodoscopes [19] measured 
two independent dE/dx and the time-of-flight (TOF) of particles. In addi- 
tion, dE/dx in the drift chamber gas was obtained as a truncated mean of 
the integrated charges of the hit-pulses. The scintillator hodoscopes consisted 
of ten upper and twelve lower plastic scintillators and photo-multiplier tubes 
(PMT's) attached at each end of scintillators. The hodoscopes were placed 
at outer-most radii and the timing resolution of each hodoscope was 55 psec 
rms, resulting in resolution of 0.014, where (3 is defined as particle veloc- 
ity divided by the speed of light. Furthermore a Cherenkov counter [20] with 
silica- aerogel radiator [21] was installed below the upper TOF hodoscopes. We 
selected the radiator having a refractive index of 1.020 (while 1.032 for BESS 
1997), in order to veto e~ / fi~ backgrounds up to 4.2 GeV. 

In 1998 flight, the payload was launched on July 29 from Lynn Lake, Manitoba 
to Peace River, Albata in northern Canada, where the cut-off rigidity is less 
then 0.5 GV/c. The scientific data were taken for a live time of 60462 sec at 
altitudes ranging from 37.4 to 35.7 km (with residual air of 4.8-5.3 g/cm 2 ). All 
the detector components worked well during the flight and 19 million cosmic- 
ray events were recorded onto two magnetic tapes (the total recorded data 
size was 33.5 Gbytes). The first-level trigger was provided by a coincidence 
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between the top and the bottom scintillators, with the threshold set at 1/3 
of the pulse height from minimum ionizing particles. The second-level trigger, 
which utilized the hit-patterns of the scintillator-hodoscopes and the inner 
drift chambers (IDC), first rejected unambiguous null- and multi-track events 
and then made a rough rigidity-determination to select negatively-charged 
particles predominantly. In addition, one of every 60 first-level triggered events 
were recorded in order to build a sample of unbiased triggers. 



3 Analysis 

Since the BESS instrument had highly symmetrical detector-configuration, we 
assumed that p's behaved similarly as protons in the instrument except for 
deflection in the magnetic field and the inelastic interactions. Thus, all se- 
lection criteria were defined based on the measured properties of protons. At 
first, we selected events with (i) a single down-going particle fully contained in 
the fiducial region of the tracking volume, (ii) only one hit in the upper TOF 
hodoscope and (iii) one or two hits in the lower TOF hodoscope. Then the 
following cuts were applied in order to ensure correct measurements. (1) The 
numbers of used hits in the refitting and in the ^-fitting should be more than 
10 and 6, respectively: (2) The reduced \ 2 °f the fitted r<p and z track had 
to be less than 6.5 and 6.0, respectively: (3) Since signals from some PMT's 
were somewhat noisy, we rejected wrong hits in the TOF hodoscopes for bet- 
ter TOF measurement: (4) The z position determined from the left-right time 
difference measured by the PMT's should match the z-impact point of the 
extrapolated track at the TOF hodoscopes: (5) The ratio of the signal ampli- 
tude of the left and right PMTs should be consistent with the z-impact point 
of the extrapolated track: The instrument was suspended with the balloon 
via a parachute (50m H x lm 0) of 340 kg, situated 24 meter above the in- 
strument. Particles passing through the parachute would interact and change 
their energy. Therefore, (6) the extrapolated track should not pass through the 
parachute: (7) The three dE/dx measurements in the scintillator-hodoscopes 
and JET chamber were loosely required as functions of R to be compatible 
with proton or p: (8) The proper value of 1/(3 was required as a function of R: 
The combined efficiency of these off-line selections was 71 - 75 % for R from 
0.5 to 5 GV/c. 

These simple and highly-efficient selections were sufficient for a very clean 
detection of p's in the low velocity ((3 < 1/1.15) region. At higher velocities, 
the e~ I n~~ background started to contaminate the p band, where we required 
the Cherenkov veto; i.e., (1) the particle trajectory to traverse the fiducial 
volume of the aerogel, and (2) the Cherenkov output to be less than 0.09 
of the mean output from e~. This cut reduced the acceptance by 20 %, but 
rejected e~ / 'fx~ backgrounds by a factor of 4400, while keeping 92 % efficiency 
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for protons and p's which crossed the aerogel with rigidity below the threshold 
(4.7 GV/c). 

Fig. 2 shows the (3' 1 versus R plot for the surviving events. We see a clean 
narrow band of 384 p's at the exact mirror position of the protons. The p band 
was slightly contaminated with the e~ / /i~ backgrounds due to the inefficiency 
of the aerogel Cherenkov counter, and in the flux calculation we subtracted 
this background whose amount was estimated as follows: (1) We counted the 
number of e~ / n~ events overlapping with the p band before the requirement 
on the Cherenkov output. (2) We obtained the probability that e~/fi~ gave 
a lower Cherenkov output than the threshold, from the distribution of the 
output for high rigidity protons (> 25 GV) which should emit enough output. 
(3) From the above two quantities we estimated the amount of the e~ / fi~ 
background as %, 0.6 % and 4.2 % at 0.25, 2 and 4 GeV, respectively. 
Backgrounds of albedo and of mis-measured positive-rigidity particles were 
totally excluded by the excellent and R~ l resolutions. To check against 
the "re-entrant albedo" background, we confirmed that the trajectories of all 
p's can be traced numerically through the Earth's geomagnetic field back to 
the outside of the geomagnetic sphere. 

We obtained the p energy spectrum at the top of the atmosphere (TOA) in the 
following way: The geometrical acceptance of the spectrometer was calculated 
by a Monte Carlo simulation [22], which was consistent with an analytical cal- 
culation [23]. The fraction of the live data-taking time was directly measured 
as 86.4 % by counting 1 MHz clock. The efficiencies of trigger and of the off-line 
selections were determined by using protons in the unbiased trigger samples. 
The TOA energy of each event was calculated by tracing back the particle 
through the detector material and the air. The survival probability of the p's 
through the air and instrument was evaluated by GEANT/GHEISHA simula- 
tion [22], which incorporated detailed material distribution, realistic detector 
performance and correct p- nuclei cross sections [24]. The e~/fi~ background 
was subtracted, whose amount is described above. And we subtracted the ex- 
pected number of atmospheric p's, produced by the collisions of cosmic rays 
in the air. The subtraction amounted to 17 ± 3 %, 21 ± 4 % and 22 ± 6 % at 
0.3, 0.7 and 2 GeV, respectively, where the errors correspond to the maximum 
difference among the three independent recent calculations [25-27]. Table 1 
contains resultant 1998 p fluxes and p/p flux ratios at TOA. The first and the 
second error represent the statistical [28] and systematic errors, respectively. 
The dominant systematic error at low energies is the uncertainty in the eval- 
uation of the interaction losses, to which we attribute ± 15 % relative error. 
At high energies, the uncertainty in the atmospheric p calculation becomes 
dominant. The statistical errors are always larger than the systematic errors 
by a large factor especially at the low energies. 
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4 Result and Conclusion 



Fig. 3 shows the BESS 1998 p spectrum, together with the previous measure- 
ments [1,10-14]. By using 1998 flight data, we have measured the p spectrum 
in a wider energy region of 0.18 to 4.20 GeV and detected again a clear peak 
around 2 GeV as we did in 1997. As recently pointed out [5], increase of solar 
activity suppresses any soft primary p component while only modestly affects 
the secondary p spectrum. The solar activity at the time of the BESS 1998 
flight was greater than the 1995 and 1997 flights [29]. Therefore, the shape of 
the 1998 spectrum should be closer to the shape of the spectrum composed 
only of secondary p, even if the primary p's existed. Shown also in Fig. 3 is 
a recent theoretical calculation [30] for the secondary p at the BESS 1998 
flight (solar modulation parameter 98 =61OMV, which was determined from 
BESS 1998 proton spectrum [33]). In the whole energy region, this calculation 
agrees well with our spectrum at 86% confidence level. This implies that sec- 
ondary p is the dominant component of the cosmic-ray p's and the model used 
for the calculation is basically correct. The expected spectrum at BESS 1997 
was also calculated in the same way by using the solar modulation parameter 
97 =5OOMV as represented by the dashed curve in Fig.3. It appears to be 
insensitive to the change of the modulation parameter. This is because the 
interstellar p flux steeply drops toward the low energy in contrast to the inter- 
stellar proton spectrum [34]. In fact, the measured spectra are nearly constant 
in the peak region. However, in the low energy region (0.18 - 0.78 GeV), the 
agreement between the 1997 spectrum and the calculation is less consistent. 
Moreover, the 1995 spectrum shows larger deviation from the calculation in 
spite of similar modulation parameter (095=54OMV) as the BESS 1997. This 
might be due to statistical fluctuation, or might suggest a contribution of low- 
energy primary p from novel source such as annihilation of neutralinos [2,4,5] 
or the evaporating primordial black holes [6,7]. Because p's from these "pri- 
mary" sources, if they exist, are expected to be prominent at low energies and 
to exhibit large solar modulations [5]. 

In order to study the solar modulation, the p/p ratio is more adequate, which 
is expected to clarify the charge-dependence of the solar modulation [3]. In 
Fig. 4 the BESS 1998 p/p ratio shows almost similar spectrum as the spec- 
trum at the solar minimum (the 1997 data), especially around 2 GeV where 
we have better statistics. The lowest energy point of our data might suffer 
some uncertainties in the proton flux due to large subtraction of atmospheric 
secondary protons, which is significant in the low energy region [33]. Despite 
the general advantage of taking p/p ratio to cancel various systematic errors, 
such as live-time and geometrical acceptance, it should be noted that the most 
significant systematic errors in the analysis come from the uncertainties in the 
atmospheric p calculation and in the interaction losses, which can not be can- 
celed even in the form of the ratio. Taking the charge- dependence of the solar 



6 



modulation into account [3], (i) the p/p ratio is expected to be nearly identical 
during the positive Sun's polarity phase; and (ii) when the polarity switches 
into the negative, the p/p ratio shall rapidly increase as represented by the 
dotted curve in Fig. 4. The BESS 1993, 1995, 1997 and 1998 spectra, measured 
in the positive polarity phase, were nearly identical as shown in Fig. 4 and 
consistent with the prediction. It can be more clearly understood in the form 
of the ratios of p/p ratio as shown in Fig. 5. The ratios of p/p ratio normalized 
by the BESS 1997 data are consistently distributed around the unity within 
the statistical fluctuation. 

Our measurements on the p/p ratio are consistent with the prediction taking 
the charge dependent solar modulation into account. However, it is too early 
to conclude that the evidence for the charge dependence is exhibited. We are 
planning further annual flights in the following solar-maximum period in which 
the polarity switches into the negative. The observation of the increase would 
be a definite evidence for the importance of the charge dependent effect in the 
solar modulation. At the same time, measurements of the p spectrum through 
the solar-maximum period itself have a crucial importance to determine the 
origin of the low energy p's. As previously mentioned, any soft p component 
from the primary sources would be suppressed as the solar activity increases, 
while the secondary p spectrum is affected modestly. Therefore, to observe 
variations in the low energy region would help us to draw a firm conclusion 
on whether we are seeing the primary p's from novel sources. 
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Table 1 

Antiproton fluxes (in xl0~ 2 m _2 s _1 sr _1 GeV _1 ) and p/p ratios (in xl0~ 5 ) at 
TOA. T (in GeV) define the the kinetic energy bins. Np and Tp, respectively, are 
the number of observed antiprotons and their mean kinetic energy in each bin. 



T (GeV) 



0.18 
0.28 
0.40 
0.56 
0.78 
0.92 
1.08 
1.28 
1.52 
1.80 
2.12 
2.52 
3.00 
3.56 



0.28 
0.40 
0.56 
0.78 
0.92 
1.08 
1.28 
1.52 
1.80 
2.12 
2.52 
3.00 
3.56 
4.20 



BESS 1998 



p flux 



p/p ratio 



2 0.24 

6 0.36 

15 0.52 

22 0.68 

23 0.84 
25 1.00 
21 1.20 

41 1.42 
29 1.66 
37 1.94 

42 2.28 
46 2.80 
46 3.30 
29 3.83 



0.34 



+0.48+0.06 
-0.23 -0.06 



n 7 ,,+0.48 +0.09 
u - '^-0.32 -0.09 



1.30 



+0.44+0.14 
-0.37-0.14 



1.23 



+0.38+0.13 
-0.33 -0.13 



9 1Q+0.60+0.21 
z - iy -0.53-0.21 



o o n +0.65+0.21 
z - JU -0.52 -0.21 



1.76 



+0.58+0.16 
-0.52 -0.16 



o on+0.65+0.29 
°- ou -0.58 -0.29 



1 oi +0.49 +0.20 
1 - Oi -0.43 -0.20 



9 9O+0.54+0.23 
Z - ZO -0.46 -0.23 



9 r.7+0.50+0.18 
z - u ' -0.36 -0.18 



9 nn +0.45+0.22 
z - uu -0.33 -0.22 



1 oq+0.45 +0.25 
1109 -0.29 -0.25 



0.20 



+0.29+0.04 
-0.14-0.04 



n .9+O.27 +0.06 
u -^ z -0.18-0.06 



n o fi +0.29+0.13 



0.95 



+0.29+0.14 
-0.26-0.14 



-, Q 7 +0.54+0.27 
i - y ' -0.48-0.27 



2.38 



+0.67+0.32 
-0.54-0.32 



2.25 



+0.74+0.31 
-0.66-0.31 



a nq+0.97+0.65 
^• 9O -0.87-0.65 



o nn+0.84+0.46 
°- uy -0.73-0.46 



k 99+1.26+0.74 
d - zz -1.08-0.74 



cr 74 +1.38+0.77 
°-' ^-1.00 -0.77 



7.201} 



+1.63+1.06 
19-1.06 



q 94+2.21+1. 

9 - z4 -1.44-l 



54 
54 



1.84 



+0.58+0.26 
-0.32 -0.26 



11 cq+3.66+2.01 
ll.DO_ 2 .05-2.01 



10 



Aerogel 



JET drift 
chambe 



Coil 




Pressure 
Vessel 



IDC 



TOF hodoscope 



1.7 m- 



Fig. 1. Cross-sectional view of the BESS 1998 detector with one of the p events. 
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Kinetic Energy (GeV) 



Fig. 3. The BESS 1998 antiproton spectrum at the top of the atmosphere, together 
with the previous data. The thick solid curve represents the expected spectrum 
for secondary p [30] at the 1998 flight (0 98 =61OMV) and well reproduces the BESS 
1998 spectrum in the whole energy region. The expected spectrum at the 1997 flight 
(</>97=500MV) was also calculated in the same way as the 1998, represented by the 
thick dashed curve. Although the curve agrees with the BESS 1997 spectrum in the 
peak region, the agreement is less consistent in the low energy region. Also shown 
are other calculations [2,3] for secondary p at solar minimum (corresponding to the 
solar activity level at the BESS 1997 flight). 



13 



-3 

10 — 



cs 
u 

a 

i a 



Solar minimum in + phase 
Solar maximum in + phase 
Solar maximum in - phase 




■ BESS(98) 
O BESS(97) 
A BESS(95) 



BESS(93) 

□ CAPRICE 

•k IMAX 

X PBAR 

LEAP 



* MASS2 



Kinetic Energy (GeV) 



Fig. 4. Comparison of the BESS 1993, 1995, 1997 and 1998 p/p ratios with previous 
data [12], and the calculation [3] taking the charge dependence of the solar modula- 
tion into account. The solid and dashed curves represent the expected p/p ratio at 
the solar minimum and at the solar maximum in the positive Sun's polarity. When 
the polarity switches, the ratio is expected to increase as represented by the dotted 
curve. 
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Fig. 5. The ratios of the p/p ratio normalized by the BESS 1997 data are distributed 
around the unity within the statistical fluctuation. The BESS 1993, 1995, 1997 and 
1998 data were measured in the positive polarity phase. Therefore, it is shown that 
the p/p ratio is nearly identical during the positive polarity phase. 
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